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phase  equations  until  a  steady  state  solution  Is  reached. 
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gas  momentum  exchange  field.  Solution  of  the  gas  phase  equa¬ 
tions  by  the  simplified  MAC  (SMAC)  method  then  completed  a 
single  cycle  of  the  iterative  procedure. 


Calculations  of  the  flow  field  due  to  an  axisymmetric 
spray  shows  good  agreement  with  experiment.  The  degrading 
effect  of  wind  on  a  spray  curtain  Is  demonstrated,  though  calcu¬ 
lations  show  that  curtain  performance  may  be  Improved  by  tilting 
the  nozzle  Into  the  wind.  The  flow  Inside  a  cooling  tower  and 
the  resulting  pressure  losses  are  seen  to  be  reasonably  predicted. 
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A  two  dimensional  model  (axi symmetric  or  planar)  has 
been  implemented  to  study  the  gas-droplet  flow  due  to  water 
sprays  or  in  cooling  towers. 

The  model  consists  of  coupled  sets  of  equations 
governing  the  gas  and  liquid  phases.  Modelling  the  gas  phase 
as  a  continuum  allows  use  of  the  Navier-Stokes  equations  with 
momentum  source  terms  included.  The  liquid  phase,  represented 
by  a  finite  number  of  droplet  trajectories  is  governed  by  the 
equations  of  motion  of  falling  spherical  droplets  subject  to 
gravitational  and  drag  forces. \  The  particle-gas  momentum 
exchange  and  the  influence  of  ^e  gas  flow  on  the  droplet  trajec¬ 
tories  results  in  a  coupl  irm...-df  both  of  these  sets  of  equations. 


^^Numerical  solution  of  these  equations  is  made  in  an 
iterative  fashion  by  solving  first  the  liquid  phase  equations 
and  then  the  gas  phase  equations  until  a  steady  state  solution 
is  reached.  Solution  of  the  liquid  phase  equations  by  a  fourth 
order  Runge-Kutta  method,  allows  subsequent  determination  of  the 
droplet-gas  momentum  exchange  field. >.  Solution  of  the  gas  phase 
equations  by  the  simplified  MAC  (SMACn  method  then  completed  a 
single  cycle  of  the  iterative  procedi^. 

^Calculations  of  the  flow  field  due  to  an  axisymmetric 
spray  shows  good  agreement  with  experiment.  The  degrading 
effect  of  wind  on  a  spray  curtain  is  demonstrated,  though  calcu¬ 
lations  show  that  curtain  performance  may  be  Improved  by  tilting 
the  nozzle  into  the  wind.  The  flow  inside  a  cooling  tower  and 
the  resulting  pressure  lo^es  are  seen  to  be  reasonably  predicted 
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1.  INTRODUCTION 

Liquid  sprays  and  cooling  towers  represent  just  a  few 
typical  examples  of  applications  of  gas-droplet  flows.  Optimi¬ 
zing  the  performance  in  these  applications  requires  that  the 
fluid  dynamics  be  understood.  Difficulties  in  modelling  such 
flows  experimentally  at  below  full-scale  levels  and  the  expense 
of  full  scale  tests  has  mandated  a  need  for  the  development  of 
numerical  models. 

The  presence  of  droplets  in  these  flow  fields  intro¬ 
duces  additional  complications  in  their  physical  descriptions 
and  subsequent  calculation.  Often  the  primary  effect  of  the 
droplets,  due  to  the  aerodynamic  drag  they  experience,  is  to 
transfer  momentum  to  the  gas,  a  consequence  of  the  conservation 
of  momentum.  This  results  in  a  modification  of  the  gas  flow 
field.  The  complications  arise  in  that  the  droplet  motion  and 
resulting  momentum  exchange  is  a  function  gas  motion  which  the 
droplets  influence.  Thus  the  complete  description  of  the  phe¬ 
nomenon  involves  coupling  between  the  two  phases. 

The  droplets  may  also  modify  other  properties  of  the 
gas,  such  as  its  turbulence  level,  density  and  viscosity,  though 
these  effects  are  treated  only  very  superficially  in  this 
report . 


Modelling  of  gas-droplet  flows  has  been  made  in  the 
past.  A  one  dimensional  axisyrametric  spray  model  has  been 
developed  at  VKI  (Ref.  13)  and  has  shown  good  results,  though 
its  applicability  for  the  inclusion  of  boundary  effects  is 
limited.  Two  dimensional  spray  and  cooling  tower  models  exist¬ 
ing  in  the  literature  (Refs.  1,2,3,8,16)  are  better  able  to 
handle  the  effects  of  boundaries.  In  this  light,  it  is  of 
interest  to  develop  a  two  dimensional  model  at  the  VKI  to  com¬ 
plement  the  range  of  applicability  of  the  one  dimensional 
axisymmetric  spray  model,  to  examine  a  planar  representation  of 
spray  curtain,  and  to  investigate  flow  in  a  cooling  tower. 
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In  this  report  the  model  used  is  outlined,  the 
numerical  method  used  for  its  solution  is  described,  and 
finally  the  results  for  the  three  applications  mentioned  above 
are  presented  and  discussed. 
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2.  THE  GAS-DROPLET  MODEL 

2 . 1  Introduction 

The  gas-droplet  model  consists  of  two  distinct  sets- 
of  equations,  one  set  governing  the  gas  phase  and  another 
governing  the  liquid  phase.  Linkage  between  these  two  sets  of 
equations  accounts  for  the  following  mode  of  gas-droplet  inter¬ 
action.  The  viscous  nature  of  the  gas  causes  droplets  moving 
relative  to  the  gas  to  experience  an  aerodynamic  drag.  Because 
momentum  must  be  conserved,  this  aerodynamic  drag,  acting  to 
change  the  momentum  of  the  droplets,  also  changes  the  momentum 
of  the  gas  in  a  signwise  opposite  manner.  The  result  is  that 
the  motion  of  the  droplets  is  influenced  by  the  local  gas 
velocity  and  visa  versa.  It  should  be  noted  that  the  effect  of 
turbulence  on  the  droplets  is  not  taken  into  account  and  the 
turbulence  induced  by  the  droplets  or  by  other  sources  is  only 
crudely  accounted  for. 

2.2  The  gas  phase  model 

The  gas  phase,  occupying  the  most  significant  portion 
of  the  flow,  is  treated  as  a  continuum.  Making  the  standard 
assumptions  of  incompressible,  isothermal,  Newtonian  fluid,  and 
assuming  the  volume  occupied  by  the  droplets  is  negligible 
allows  the  following  form  of  the  Navier-Stokes  equations  to  be 
applied.  The  possibility  of  considering  planar  or  axisymmetric 
cases  is  provided. 
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Here  V^,  V^,  P  are  the  horizontal  and  vertical  components  of 
velocity  and  pressure  respectively,  and  r,  z  are  the  horizontal 
and  vertical  coordinates,  p  and  v^pp  represent  the  gas  density, 
and  the  effective  viscosity  of  the  gas  which  could  crudely 
account  for  the  turbulence  induced  by  the  droplets,  a  is  the 
switch  for  the  axisymmetric  (a=l)  or  planar  (a=0)  cases. 


F  and  F  are  the  horizontal  and  vertical  momentum 
r  z 

source  terms  which  represent  external  sources  (sinks)  of  gas 
momentum  per  unit  volume  such  as  that  due  to  the  droplets.  They 
are  defined  as  the  sum  of  the  external  contributions  to  the 
horizontal  or  vertical  components  of  gas  momentum  inside  a 
volume  aV  divided  by  the  volume  aV ,  expressed  symbolically 
bel ow. 
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The  horizontal  and  vertical  momentum  source  terms  due 
to  the  presence  of  droplets  may  be  obtained  by  summing  respec¬ 
tively  the  horizontal  and  vertical  components  of  aerodynamic 
drag  acting  at  each  droplet  inside  a  volume  aV  and  dividing  by 
the  volume  aV. 


These  equations  and  appropriate  boundary  conditions 
are  applied  inside  a  representative  domain  (see  Fig,  2.1). 

Four  sets  of  boundary  conditions  are  considered  here 
corresponding  to  the  following  types  of  boundaries no-sl i p 
walls,  slip  walls,  specified  inflow/outflow  boundaries,  and 
unspecified  inflow/outflow  boundaries. 


A  no-slip  boundary  condition  is  applied  in  the  presence 
of  a  rigid  wall  where  the  friction  effects  of  such  a  boundary 
are  important.  The  component  of  velocity  tangent  to  the  wall  is 
set  equal  to  zero  satisfying  the  no-slip  condition  while  the 
normal  component  is  set  to  zero  to  satisfy  the  condition  of  zero 
mass  flow  across  the  rigid  wall.  The  pressure  boundary  condi¬ 
tion  follows  by  applying  the  momentum  equation  in  the  direction 
normal  to  the  wall.  These  conditions  are  expressed  below.  Here 
the  subscript  N  denotes  the  normal  coordinate  while  t  denotes 
the  tangential  coordinate. 
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The  slip  wall  boundary  condition  is  applied  at  an 
axis  of  symmetry  or  at  a  rigid  wall  where  the  friction  effects 
of  such  a  boundary  are  not  important  or  not  desired.  The  compo¬ 
nent  of  velocity  normal  to  the  wall  is  set  equal  to  zero  to 
again  satisfy  the  condition  of  zero  mass  flow  across  the  rigid 
wall.  The  slip  condition  is  satisfied  by  requiring  the  gradient 
of  the  tangential  component  to  be  zero  at  the  wall.  Again  the 
pressure  boundary  condition  follows  from  the  momentum  equation 
in  the  direction  normal  to  the  wall. 


V  =0 
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— t  =  0  (2.7) 
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The  boundary  conditions  for  the  velocity  at  a  specified 
inflow/outflow  boundary  follows  by  definition,  while  the  pressure 
boundary  conditions  may  be  obtained  by  substituting  the  specified 
velocities  into  the  momentum  equation  in  the  direction  normal  to 
the  boundary. 


*  V^(r.z) 
=  V^(r,z) 
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—  -►  obtained  from  momentum  equation, 
a  N 

At  the  unspecified  inflow/outflow  boundary  the  applied 
boundary  conditions  reflect  the  porous  nature  of  such  a 
boundary . 


=0  (2.9) 

■*  obtained  from  momentum  equation 

It  should  be  noted  that  since  the  pressure  boundary 
conditions  are  a11  von  Neumann  type,  the  pressure  needs  to  be 
specified  at  one  point  In  the  domain. 

2.3  The  droplet  phase  model 

The  droplet  phase  Is  modelled  by  considering  a  finite 
and  discrete  distribution  of  droplets  of  varying  diameter, 
injection  velocity,  and  injection  angle.  Using  a  Lagrangian 
approach  Individual  droplets  are  followed  from  injection  until 
hitting  the  ground.  Appropriately  assigning  a  portion  of  the 
total  number  of  droplets  Injected  by  nozzle  per  unit  time  to 
each  of  the  trajectories  and  determining  the  drag  force  along 
each  trajectory  allows  representation  of  the  droplet-gas 
momentum  exchange.  This  approach  assumes  no  droplet-droplet 
1  nteractioi'iS  such  as  collisions  or  droplet  break-up  occur. 

Thus  the  equations  of  motion  of  the  droplet  need  to 
be  considered.  Forces  acting  on  the  droplet  other  than  aero¬ 
dynamic  drag  and  gravitational  force  are  neglected.  Further 
assuming  that  the  droplets  are  spherical  and  non-evaporating 
allows  the  following  equations  of  motion  to  be  written  : 
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du 

m  -  »  -  f  ±  mg  (2.11) 


dt 


Here  u^.,  are  the  horizontal  and  vertical  components  of  drop¬ 
let  velocity  and  R,Z  are  the  horizontal  and  vertical  coordinates 
of  the  droplet  trajectory.  The  mass  of  the  drop!  et ,  ‘  "ni" ,  is 
simply  defined. 


m  =  p  - 
«  6 


(2.14) 


where  D  is  the  droplet  diameter  and  is  the  droplet  density. 

The  horizontal  and  vertical  components  of  aerodynamic 

drag  force  f_,  f,  needed  for  the  calculation  of  and  F,  may 
r  2  r  2 

be  written  as  follows 


f,  =  C,Re 


(2.15) 


^  =  CjjRe  —  (u^-v^) 


(2.16) 


The  Reynolds  number  of  the  droplet.  Re,  is  defined. 


(Uz-V2)^+(u^-v^)^  D 


(2.17) 


The  drag  coefficient  of  the  droplet  is  assumed  equiva' 
lent  to  that  of  a  sphere  and  is  described  by  the  following 
standard  form  fit  (Ref.  9) 
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Cn  =  —  +  +  .4  0.  <  Re  <  10®  (2.18) 

^  Re  iWRe 

Initial  conditions  : 


The  initial  conditions  for  the  droplet  injection  are 
obtained  by  considering  the  initial  position,  injection  velo¬ 
city,  and  injection  angle  (see  Fig.  2.1).  These  variables  are 
dependent  on  the  particular  phenomenon  being  studied  such  as 
sprays  where  the  droplets  are  injected  by  a  nozzle  or  cooling 
towers  where  the  droplets  drip  from  the  packing. 


R(t=0) 
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U^(t=0) 

=  UoSine 

U^(t^O) 
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Here  R^ 

is  the  horizontal 

distance  from  the  left  boundary  of  the 

domain  to  the  injection  point,  is  the  vertical  domain  length, 
H  is  the  height  from  the  ground  to  the  injection  point,  Uq  is 
the  injection  velocity  and  0  is  the  injection  angle. 
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3.  THE  NUMERICAL  METHOD 

It  is  difficult  indeed  to  conceptualize  a  numerical 
method  for  the  simultaneous  solution  of  both  sets  of  equations. 
A  more  natural  approach  involves  a  two-step  iterative  procedure 
where,  in  the  first  step,  the  solution  of  the  gas  phase  equa¬ 
tions  is  made  using  the  most  updated  solution  of  the  droplet 
equations  and,  in  the  second  step,  the  solution  of  the  droplet 
equations  is  made  using  the  newest  gas  phase  solution.  This 
process  continues  until  convergence  of  both  sets  of  equations 
is  reached. 


Such  a  procedure,  common  in  the  literature  for  gas- 
droplet  flows  (Refs.  1,2,3),  is  used  here  for  the  solution  of 
the  gas-droplet  model  equations.  The  numerical  methods  used  to 
solve  each  set  of  equations  reflect  the  different  nature  of  the 
equations  themselves  and  are  discussed  separately  in  the  follow¬ 
ing  two  sections.  The  third  section  describes  the  calculation 
of  the  source  terms  and  finally  in  a  fourth  section  the  details 
of  the  two-step  procedure  used  here  are  discussed. 

3 • 1  Numerical  solution  of  the  gas  phase  equations 

The  numerical  scheme  for  the  solution  of  the  gas 
phase  equations  is  based  on  the  Simplified  Marker  and  Cell 
(SMAC)  method  (Refs.  4,5).  The  SMAC  method,  used  for  the 
solution  of  the  primitive  variable  (u,v,P)  form  of  the  Navier- 
Stokes  equations,  is  a  finite-difference  time-marching  procedure 
characterized  by  a  cell  structure  shown  in  figure  3.1. 

Pressure  is  defined  at  the  center  of  each  cell, 
horizontal  components  of  velocity  defined  at  the  vertical  cell 
boundaries,  and  vertical  components  of  velocity  defined  at  the 
horizontal  boundaries  of  the  cell.  The  vertical  and  horizontal 
momentum  source  cells  are  rectangular,  of  dimension  Arx&z,  and 
centered  at  the  positions  where  the  corresponding  vertical  and 
horizontal  velocity  components  are  defined. 


The  method  is  an  iterative  procedure  involving  four 
steps  at  each  cycle.  During  the  first  step  the  velocity  compo¬ 
nents  inside  the  domain  are  explicitly  updated  using  a  discre¬ 
tized  form  of  the  momentum  equations.  Using  the  values  of 
velocity  inside  the  domain,  the  unspecified  normal  velocities 
on  the  domain  boundary  are  calculated  to  complete  the  second 
step.  In  the  third  step  sweeps  are  made  through  the  grid  check¬ 
ing  to  make  sure  that  the  divergence  (mass  conservation)  of 
each  cell  is  below  some  tolerable  level.  If  the  divergence  is 
greater  than  this  value  then  the  pressure  inside  the  cell  is 
modified  appropriately  and  the  velocities  are  redistributed  so 
that  continuity  is  satisfied.  The  sweeps  through  the  grid  con¬ 
tinue  until  the  divergence  of  each  cell  is  below  the  tolerence 
level  . 

Finally  in  the  fourth  step  the  boundary  conditions 
necessary  to  specify  the  velocity  components  tangent  to  the 
boundaries  in  the  "fictitious"  cells  surrounding  domain  are 
applied . 

This  four-step  procedure  is  applied  iteratively  until 
convergence  is  obtained. 

3.1.1  Eq^ation^  for  the  first  £tep 

As  mentioned  the  equations  used  in  the  first  step  are 
finite  difference  forms  of  the  momentum  equations  (2.2,  2.3). 
Here  forward  differencing  in  time,  upwind  differencing  on  the 
advection  terms  and  central  differencing  of  the  other  spatial 
derivatives  is  used  to  obtain  explicit  equations  for  the  velo¬ 
city  coiiponents  inside  the  domain.  The  form  of  these  equations 
is  similar  to  that  used  by  Viecelli  (Ref.  5),  and  are  given  as 
follows.  Here  the  tilde  (~)  represents  the  updated  velocities. 
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Axial  momentum  equation 
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Discretized  velocities  not  explicitly  available  are 
obtained  by  a  simple  averaging  such  as 
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A  flagging  scheme  is  used  to  allow  the  possibility 
of  obstacles  inside  the  domain.  The  flags  f  .  j  +  1/2 
^i+1/2  j  equations  (3.1,  3.2)  are  set  equal  to  zero  inside  or 
on  the  boundary  of  an  obstacle  and  equal  to  one  elsewhere  in 
the  flow  domain.  This  causes  the  components  of  velocity  on  or 
inside  obstacles  to  be  unchanged  during  the  first  step,  while 
allowing  the  other  velocity  components  inside  the  domain  to  be 
updated  according  to  the  discretized  momentum  equations. 


3.1.2  Equations  for  the  second  step 


After  completing  the  first  step,  the  velocity  compo¬ 
nents  inside  the  flow  domain  will  have  been  updated.  In  the 
second  step,  the  unspecified  normal  velocity  components  on  the 
boundary  are  updated  by  requiring  that  the  volume  flow  through 
the  boundary  be  the  same  as  that  on  the  cell  face  just  inside 
the  domain.  Normal  velocity  components  on  boundaries  which 
have  specified  inflow/outflow,  such  as  non-porous  walls, 
symmetry  boundaries,  or  boundaries  where  the  inflow/outflow  is 
known,  are  left  unchanged  by  this  step.  Again  making  use  of 
the  flagging  scheme  the  equations  of  the  second  step  may  be 
written.  Here  only  the  equations  for  the  top  and  right  boun¬ 
daries  are  given  because  extension  to  the  bottom  and  left  boun¬ 
daries  is  straightforward. 
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where  NX  is  the  number  of  cells  in  the  horizontal  direction 
including  the  two  fictitious  cells. 


3.1.3  2he  equat^ions  for  th£  thi^r^ 

Arriving  at  the  third  step,  all  the  velocity  compo¬ 
nents  inside  the  domain  and  on  the  domain  boundaries  have  been 
updated.  Now  the  divergence  of  each  cell  is  checked  and  the 
pressure  modified  and  velocities  reapportioned  if  the  divergence 
is  greater  than  a  specified  value. 
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yields  the  following. 
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Pres^ure_and_vel oci^t^  modif i cation 


If  D..  is  greater  than  a  specified  tolerance  then  the 
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following  operations,  similar  to  that  used  by  Hirt  and  Cook 
(Ref.  4),  are  applied  to  update  the  pressure  and  velocity  in 
the  cell.  Updated  quantities  are  devoted  by  a  tile  (~). 
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e  is  the  relaxation  coefficient  and  varies  between  0.  and  2. 
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The  process  of  checking  the  divergence  continues  cell 
by  cell  through  the  mesh  until  a  complete  sweep  is  made  without 
any  cell  having  a  divergence  greater  than  the  tolerence. 


This  tolerance,  e,  is  typically  .0001-. 001  ^  where 
U  is  a  characteristic  velocity  and  L  a  length  scale  of  the  flow. 

3.1.4  Equations  for  the  fourth_ste£ 

The  fourth  step  consists  of  specifying  the  velocity 
components  tangent  to  the  boundary  in  the  ficittious  cells 
which  surround  the  domain  and  inside  the  cells  which  represent 
obstacles  in  order  to  satisfy  the  remaining  boundary  conditions. 
Two  types  of  boundary  conditions  for  velocity  components  tangent 
to  the  boundaries  are  considered;  those  involving  slip  walls 
and  those  involving  no-slip  walls. 

No^sii£  wal^l^ 

The  no-slip  wall  is  a  rigid  wall  where  there  is  no 
tangential  component  of  velocity  at  the  wall  and  is  often  used 
to  represent  ceilings,  ground  etc.,  where  the  friction  effects 
of  such  boundaries  are  important.  Because  the  boundary  of  the 
domain  coincides  with  the  cell  boundaries  and  because  of  the 
staggered  grid  used,  averaging  of  the  velocity  components 
tangential  to  the  wall,  located  a  half  cell  width  inside  and 
outside  the  boundary,  is  required  to  satisfy  the  no-slip 
conditions 

V-l*2  *  ■  ’w  '  “ 

Vl/2  •  -  »„*l/2 

domain  boundary 
fictitious  cells 
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Slj^p^wal  l^s 


The  slip-wall  condition  may  be  applied  when  no  tan¬ 
gential  velocity  gradient  at  the  boundary  is  desired,  as  for 
the  cases  of  axis  of  symmetry  boundaries,  unspecified  inflow- 
outflow  boundaries,  or  rigid  walls  where  the  friction  effects 
are  not  of  interest. 

Again  averaging  across  the  boundary  is  used  to 
obtain  the  following. 

V  -V 

wtl/2  vv-1/2  ^  Q 
AZ 


V  =  V 

w-1  +  2  ’'w+1/2 


(3.16) 


^ojiver^ence  ter_ia 

The  test  for  convergence  of  the  gas  phase  is  made  by 
examining  the  residual  of  the  velocities  defined  below  and  by 
examining  the  velocities  and  pressure  at  two  points  in  the 
flow  field. 


N 

N-1  ; 

i. 

r  N 

N-1 

8  =  22 

V. 

+  22 

V 

i  j 

ij+1/2 

^ij+l/2j 

ij 

.  '"i  +  l/2,j 

^i  +  1  ,j. 

A  typical  value  of  the  residual  at  convergence  is  on 

-7  “2  " 

the  order  of  1.  x  lo  *  U  where  U  is  a  characteristic  velocity 
]iume£i  cal^ 

From  references  4  and  5  the  following  two  stability 
criteria  are  deduced. 

First  criteria  is  a  result  of  the  finite  difference 
approximation  of  viscous  terms  in  the  momentum  equations. 


At  < 
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(3.18) 
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The  second  criteria  requires  that  a  fluid  particle 
not  cross  a  cell  in  a  single  time  step. 


At  <  min 


(3.19) 


3 . 2  Numerical  solution  at  the  droplet  equations 

The  droplet  equations  are  solved  simultaneously  using 
a  fourth-order  Runge-Kutta  technique  for  each  different  set  of 
initial  conditions.  Such  techniques  are  well  documented  in  the 
literature  (Ref.  6)  and  thus  the  details  are  not  considered  here 

One  detail  that  needs  to  be  clarified  concerns  the 
interpolation  of  the  gas  velocities  which  are  not  always  avail¬ 
able  at  the  same  position  in  the  flow  field  as  the  droplet. 

A  bi-variate  linear  interpolation,  similar  to  that  used  in  MAC 
codes  to  calculate  the  marker  particle  positions,  is  applied 
here  (Ref,  7).  Thus,  the  marker  particles  of  the  original  MAC 
method  can  be  thought  of  as  merely  massless  "droplets"  which 
do  not  interact  with  the  gas  as  opposed  to  the  droplets  in  the 
method  discussed  here  which  do. 


Once  the  droplet  equations  have  been  solved,  the 
solution,  the  droplet  position  and  velocity  along  the  trajec¬ 
tories  as  a  function  of  time,  may  then  be  used  to  calculate 
the  droplet  to  gas  momentum  source  terms. 

3.3  Calculation  of  the  droplet  to 
gas  momentum  source  term 

Once  a  droplet  trajectory  has  been  calculated  the 
contribution  of  the  droplets  along  the  trajectory  to  the 
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droplet-gas  momentum  exchange  may  be  evaluated.  The  trajectory 
consists  of  a  number  of  segments  the  location  of  whose  end 
points  have  been  computed.  The  horizontal  and  vertical  compo¬ 
nents  of  the  droplet  velocity  are  also  known  at  these  end  points 
and  allow  the  aerodynamic  drag  which  acts  on  the  droplets  along 
the  trajectory  to  be  determined. 

The  horizontal  and  vertical  droplet  to  gas  momentum 
source  terms  may  be  determined  from  the  sum  of  the  respective 
components  of  the  aerodynamic  drag  acting  on  the  droplets  in 
the  source  cells,  divided  by  the  volume  of  the  source  cells 
as  defined  below. 


.j+1/2 


'^i  .j  +  1/2 
k 


^•j  +  1/2 


or  equivalently 


■i. j  +  1/2 


AV 


i. j+1/2 


'■i+l/2,j 


AV 


i+l/2,j 


where  fr.„„  *nd  are  the  average  horizontal  and  vertical 

oVc  TtaVc 


aerodynamic  drig  force  acting  on  the  j  j  +  i/2 

droplets  in  tl 
respectively. 


droplets  in  the  source  cells  of  volume  AV.  .  and  aV .  .  . 


The  evaluation  of  the  droplet  to  gas  momentum  source 
proceeds  segment  by  segment  along  each  trajectory  from  the 
Injection  point  until  the  droplet  reaches  the  ground.  The 
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droplet  to  gas  momentum  source  evaluation  due  to  droplets  along 
a  single  segment  begins  by  determining 

which  of  the  source  cells  the  segment  is  contained  in  (see  Fig. 
3.2).  The  segment  is  divided  into  subsegments  if  the  segment 
crosses  any  source  cell  boundaries.  The  number  of  droplets 
along  each  subsegment  is  calculated  as  follows  : 


N 


am 

^bot 


IT 

6 


D 


3 


where  at  is  the  time  step  of  the  droplet  integration,  is 

the  length  of  the  segment  and  am  is  the 

subsegment  length.  f„^  is  the  volume  flow  of  liquid  assigned 

qw 

to  the  trajectory  of  interest,  and  is  calculated  such  that 
the  distribution  of  droplet  quantity  summed  over  all  the  tra¬ 
jectories  in  the  domain  represents  the  physical  phenomenon 
being  studied. 

The  calculation  of  the  average  force  is  made  by 
simply  averaging  the  aerodynamic  drag  force  acting  on  the 
droplet  at  both  of  the  segmant  endpoints. 

With  the  simple  determination  of  the  source  cell 
volume  the  evaluation  of  the  droplet  to  gas  momentum  source 
contribution  from  this  segmant  is  complete. 

If  two  or  more  segments  are  contained  in  the  same 
source  cell  the  contributions  of  each  to  the  momentum  source 
are  simply  added. 

It  is  to  be  noted  that  the  time  step  for  the  integra¬ 
tion  of  the  droplet  equation  should  be  small  so  that  the  dis¬ 
tance  travelled  by  the  droplet  is  small  compared  with  the  cell 
sizes,  ensuring  accuracy  in  the  evaluation  of  the  momentum 
source  terms. 
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3 , 4  The  solution  procedure 

A  solution  procedure  similar  to  the  one  discussed 
in  the  beginning  of  this  section  was  adopted  and  is  outlined 
in  figure  3.3. 

The  procedure  begins  with  a  “guessed"  initial  gas 
flow  field.  The  initial  droplet  trajectories  and  droplet  to 
gas  momentum  source  are  calculated  using  this  initial  flow 
field. 


The  gas  phase  computations  are  then  made  using  the 
newly  obtained  momentum  source  terms.  Because  the  changes  in 
the  gas  flow  field  from  cycle  to  cycle  are  small,  ten  cycles 
of  the  SMAC  procedure  are  made  before  resolving  the  droplet 
equations . 


After  updating  the  droplet  trajectories  and  droplet 
to  gas  momentum  source  terms  based  on  the  most  updated  gas 
flow  field,  the  gas  flow  computations  continue  with  10  more 
cycles . 


The  process  is  repeated  until  convergence  of  the 
gas  flow  field  is  obtained. 
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4.  APPLICATIONS  AND  RESULTS 


The  method  is  applied  to  three  typical  problems  invol¬ 
ving  gas-droplet  flows;  calculation  of  the  flow  field  induced  by 
an  axisymmetric  spray,  determination  of  the  wind  effects  on  a 
spray  curtain,  and  examination  of  the  flow  in  the  lower  portion 
of  a  cool i ng  tower . 

4 . 1  Axisymmetric  spray 

The  first  application  made  is  in  the  prediction  of 
the  flow  induced  by  an  axisymmetric  spray. 


4.1.1  ModeJ^li^n^  of^aji  ^x;[symme_tr i^c_S£r^y 

The  axisymmetric  spray  is  represented  by  a  discrete 
and  finite  distribution  of  droplets  of  varying  sizes  and  initial 
injection  angles  injected  into  the  gas  from  a  single  point  along 
the  axis  of  symmetry. 


A  distribution  of  droplet  sizes  is  considered  using 
an  approach  taken  by  Alpert  and  Mathews  (Ref.  10).  They  attri¬ 
bute  Dundas  (Ref.  11)  with  having  shown  that  the  range  of  drop¬ 
let  sizes  produced  by  a  deflector  type  nozzle  is  represented 
by  the  Rassi n-Ramml er  distribution  function  given  below. 


F{D)  =  1-exp 


2' 

-{tn2)  \ 
OZ 


(4.1.1) 


where  F(D)  is  the  mass  fraction  of  droplets  having  diameters 

less  than  D  and  0  is  the  mean  droplet  size.  Five  classes  of 
m 

droplet  sizes  are  considered,  each  having  the  same  proportion 
of  the  liquid  volume  flow.  Thus,  the  representation  diameters 
of  each  class  are  .39,  .72,  1.0,  1.32  and  1.82  times  the  mean 
droplet  diameter,  calculated  using  the  simple  approach  given  by 
Heskested  et  al.  (Ref.  12)  as  follows  : 
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(4.1.2) 


where  do  is  the  inside  diameter  of  the  nozzle  and  aP^  is  the 
water  delivery  pressure.  Due  to  the  lack  of  data,  the  nozzle 
design  parameter,  Cj^,  is  taken  as  l.G. 

Each  class  of  droplet  sizes  is  further  discretized 
by  considering  a  distribution  of  initial  injection  angles  in 
order  to  represent  a  full  cone  spray.  The  initial  injection 
angles  are  chosen  such  that  a  line  perpendicular  to  the  axis 
of  symmetry  is  cut  into  the  number  of  trajectories  per  class 
minus  one  segments  of  equal  length  by  a  set  of  lines 

tangent  to  the  initial  droplet  trajectories .  This  is  expressed 
by  the  fol lowing 

-1 

fl'j.  =  tan 

where  is  the  half  angle  of  the  spray. 

It  is  desired  to  have  a  mass  flow  of  droplets  per 
unit  area  which  is  fairly  constant  throughout  the  spray  near 
the  nozzle.  The  following  formula  are  used  to  assign  the 
volume  flow  of  liquid  to  each  trajectory. 


tane 
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i  =  l, . . .  ,N 


traj 


(4.1.3) 
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(4.1.4) 
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The  one-fifth  factor  represents  the  equal  share  of  the  total 
volume  flow  of  water  assigned  to  the  five  droplet  size  classes. 
The  ratio  on  the  right  hand  side  of  these  equations  represents 
the  area  of  the  i^*’  ring  formed  by  circles  drawn  at  the  midpoint 
of  the  segments  discussed  above,  divided  by  the  total  area  of  a 
circle  drawn  at  the  endpoint  of  the  outermost  segment.  This 
ratio  gives  a  fairly  constant  mass  flow  per  unit  area  for  mode¬ 
rate  half  angles  of  the  spray. 


A  constant  droplet  injection  velocity  is  chosen  due 
to  the  lack  of  data  concerning  this  parameter,  and  is  calculated 
by  assuming  that  the  injection  velocity  of  the  droplet  is  equal 
to  the  velocity  of  the  liquid  at  the  outlet  of  the  nozzle.  Such 
an  approach  is  used  by  Buchlin  (Ref.  13)  and  Botterrud  (Ref.  11) 


(4.1.5) 


Expressing  the  condition  that  the  droplets  are  injec¬ 
ted  from  a  single  point  along  the  axis  of  symmetry  completes 
the  specification  of  the  inlet  conditions  for  the  droplet  phase. 


=  height  of  spray  (4.1.6) 

Computations  are  made  on  a  rectangular  domain  shown 
in  figure  4.1.1  using  a  20  x  20  to  20  x  40  cell  mesh.  The 
boundary  conditions  which  are  applied  are  also  shown  in 
figure  4.1.1. 

Optimization  of  the  relaxation  coefficient,  e,  for 
the  pressure-velocity  modification  step  of  the  numerical  method 
yielded  a  value  of  1.4. 

The  value  of  the  viscosity  is  taken  as  the  molecular 
viscosity  of  air. 
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4.1.2  Results 

Calculation  of  the  flow  fields  due  to  an  axi symmetric 
spray  for  several  different  nozzles  is  made. 

The  flow  fields  due  to  a  Lechler  SZl  nozzle  with  the 
characteristics  shown  in  Table  4.1.1  are  displayed  in  figures 
4. 1.2-4. 1.7.  Here  calculations  are  performed  in  a  square  domain 
using  a  20  x  20  cell  mesh  with  the  nozzle  located  on  the 
ceiling.  The  first  three  figures  display  the  flow  field  with 
the  nozzle  two  meters  above  the  ground  for  increasing  mass  flow 
through  the  nozzle,  from  6.6  to  21.1  liters  per  minute.  These 
figures  represent  typical  nozzle  operating  modes  for  use  in 
spray  curtains  or  for  fire  extingtion.  The  second  three  figures 
display  the  flow  fields  for  the  same  corresponding  mass  flows 
but  with  the  nozzle  located  0.4  meters  above  the  ground  allowing 
the  near  field  behaviour  of  the  nozzle  to  be  evaluated.  To  be 
noted  in  all  six  of  these  figures  is  the  manner  in  which  the  air 
enters  the  spray.  Outside  the  spray  the  velocity  vectors  are 
nearly  perpendicular  to  the  axis  of  symmetry  but  are  quickly 
deflected  parallel  to  this  axis  after  entering  the  spray.  The 
entrained  air  then  leaves  the  spray  in  a  thin  layer  after  impin¬ 
ging  on  the  wall.  These  results  are  observable  in  both  experi¬ 
ment  and  in  other  numerical  calculations  (Refs.  10,15). 

Also  :5een  in  these  figures  is  the  increase  in  the 
magnitudes  of  the  velocities  inside  the  spray,  or  equivalently, 
an  increase  in  the  volume  of  air  entrained  into  the  spray  with 
increasing  volume  flow  through  the  nozzle. 

The  continuous  lines  rhown  in  these  figures  represent 
the  outermost  trajectory  of  each  of  the  fivr-  droplet  size 
classes,  with  the  trajectory  closest  to  the  axis  of  symmetry 
representing  the  smallest  droplet  size  and  the  furthest,  the 
largest  droplet  size. 
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For  the  two  meter  high  spray  the  radial  distribution 
of  droplets  inside  the  spray  is  highly  dependent  on  the  droplet 
diameter  with  the  smallest  droplets  concentrated  near  the  axis 
of  symmetry  for  the  bulk  of  the  spray.  The  reason  for  the 
decrease  in  the  penetration  of  the  smaller  droplets  is  due  to 
their  higher  drag  to  weight  ratio  and  to  the  fact  that  their 
terminal  velocities  are  well  below  the  injection  velocity. 

It  should  be  noted  that  droplet  transport  due  to  the 
turbulence  within  the  spray  and  droplet  break-up  are  not  con  i- 
dered  in  the  model,  both  of  which  could  allow  smaller  droplets 
to  migrate  further  out  from  the  axis  of  symmetry. 

For  the  0.4  meter  high  spray  and  close  to  the  nozzle 
for  the  two  meter  spray  the  droplets  are  fairly  well  distributed 
radially,  particularly  for  the  larger  droplet  sizes. 

4. 1.2.1  Comparison_and_experiment 

Comparison  of  the  calculated  amount  of  air  entrained 
into  the  spray  is  obtained  by  plotting  the  entrainment  efficiency 
versus  the  nondimensional  spray  envelope  diameter  squared.  The 
entrainment  efficiency,  n.  is  defined  as  the  ratio  of  the  volume 
flow  of  air  entrained  into  the  spray,  Q, ,  to  the  volume  flow  of 

a 

liquid  through  the  nozzle,  Q^,  while  the  inverse  spray  envelope 
diameter,  D^,  squared  is  nondimensionalized  as  given  below 


r  V^(r)dr 


X  = 


28 


Typical  results  of  experimental  investigations  carried 
out  at  VKI  (Ref.  15)  are  plotted  in  figure  4.1.8,  Superimposed 
on  the  experimental  data  is  the  McQuaid  correlation  and  the  ±20% 
limit?  between  which  much  of  the  data  is  located. 

The  calculated  values  for  the  cases  discussed  above 
are  shown  in  figure  4.1.9.  Fair  to  good  agreement  is  seen, 
particularly  further  away  from  the  nozzle.  Near  the  nozzle  the 
entrainment  seems  to  be  slightly  underpredicted. 

Slightly  better  agreement  is  seen  in  figure  4.1.10. 

Here  the  results  of  the  calculations  are  displayed  for  the  same 
nozzle  characteristics  except  th.  the  half  angle  of  the  spray 
is  reduced  from  30  degrees  to  22.5  degrees.  The  30  degree  half 
angle  represents  the  value  given  by  the  manufacturer,  while  the 
22.5  degree  half  angle  was  the  value  measured  by  Buchlin  (Ref.  15). 

Similar  calculations  are  made  for  the  SZ2  nozzle  with 
characteristics  shown  in  Table  4.1.2.  Again  the  nozzle  is  loca¬ 
ted  on  the  ceiling  and  a  squared  domain  with  a  20  x  20  mesh  is 
used.  The  values  of  the  entrainment  efficiency  are  shown  plot¬ 
ted  in  figure  4,1.11.  Good  agreement  is  depicted. 

Interestingly  the  overlap  of  the  calculations  made  on 
three  different  size  domains  shows  little  variation  in  the 
efficiency  indicating  first,  that  the  boundaries  have  little 
influence  and  secondly,  that  there  is  sufficient  resolution 
near  the  nozzle  to  calculate  the  efficiency  correctly. 

Comparison  between  the  calculated  and  experimentally 
measured  axial  velocity  profiles  within  the  spray  are  shown  in 
Figs.  4.1.12-4.1.14  for  three  different  locations  in  the  spray. 

The  numerical  values  are  taken  from  the  calculations  made  on 
the  one  meter  by  one  meter  domain.  Fair  agreement  is  seen, 
though  the  velocity  profile  seems  to  be  more  underestimated 
close  to  the  nozzle  and  further  from  the  axis  of  symmetry. 
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The  slight  underestimation  of  air  entrained  into  the 
spray  noted  previously  is  consistent  with  these  results.  Pos¬ 
sible  improvements  could  be  attained  through  the  use  of  a  tur¬ 
bulence  model  which  might  further  encourage  entrainment  of  air 
into  the  spray. 

4. 1.2. 2  Eval uationof  ceiling  effect 

The  model  is  used  to  evaluate  the  effect  of  the 
ceiling  on  the  performance  of  the  spray.  Calculations  are  made 
for  the  SZl  nozzle  with  a  volume  flow  of  21.1  liters  per  minute 
located  halfway  below  the  ceiling  and  the  ground.  Rectangular 
domains  are  used,  of  dimensions  two  meters  wide  by  4  meters 
high  and  0.4  meters  wide  by  0.8  meters  high  with  a  20  x  40  cell 
mesh.  Results  are  compared  with  the  corresponding  calculations 
with  the  nozzle  located  on  the  ceiling  and  are  shown  in 
figure  4.1.15. 

Only  very  slight  improvement  in  the  efficiency  of  the 
spray  is  seen  as  the  ceiling  is  displaced  from  the  nozzle,  with 
the  largest  differences  observed  near  the  nozzle  though  these 
differences  may  be  exaggerated  by  the  lack  of  mesh  resolution 
near  the  nozzle.  It  is  clear,  hoi/ever,  that  tne  bulk  of  the 
flow  inside  the  spray  is  only  slightly  effected  by  the  presence 
of  the  ceiling.  This  is  confirmed  by  Botterrud  (Ref.  14)  and 
Weinacht  (Ref.  16). 

4. 1.2. 3  Influence  of_pressure  terms 
momentum_eguations 

The  influence  of  the  pressure  terms  in  the  momentum 
equations  may  be  examined  by  computing  the  flow  field  with  the 
pressure  terms  eliminated  from  the  momentum  equations  and  com¬ 
paring  the  results  with  the  corresponding  calculations  discussed 
previously.  The  results  of  the  computations  for  the  SZl  nozzle 
with  30  degree  half  angle  and  21.1  liter  per  minute  volume  are 
shown  in  figure  4.1.16.  The  results  show  little  difference  in 
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the  amount  of  air  entrained  into  the  spray  for  the  calculations 
with  and  without  pressure  terms. 

This  result  is  significant  in  that  it  confirms  the 
validity  of  neglecting  the  pressure  gradient  in  the  approach 
taken  by  Buchlin  (Ref.  13). 

4.1.3  Concl usions 

The  results  show  fair  to  good  prediction  of  the  flow 
inside  the  spray  demonstrating  the  adequacy  of  the  model. 

Though  entrainment  is  slightly  underestimated  improvement  may 
be  possible  with  the  addition  of  a  turbulence  model. 

The  lack  of  influence  of  the  ceiling  on  the  entrainment 
efficiency  over  the  bulk  of  the  spray  is  demonstrated  by  the 
cal cul ations . 

Results  also  indicate  that  the  pressure  terms  in  the 
momentum  equations  are  of  little  importance  over  the  majority 
of  the  flow  field. 
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d 

0 

mm 

Qw 

5,t/min 

4Pn 

KPa 

®MAX 

degrees 

4.4 

6.6 

34.3 

30 

4.4 

13.3 

176.5 

30 

4.4 

21.1 

617.8 

30 

Table  4.1.1  Properties  of  Lechler  SZl  nozzle 


do 

Qw 

®MAX 

mm 

it/mi n 

KPa 

degrees 

5.6 

36.0 

588.6 

30 

5.6 

44.5 

981. 

30 

Table  4.1.2  Properties  of  Lechler  5Z2  nozzle 
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4 . 2  Wind  effects  on  a  planar  spray  curtain 

As  a  second  application  of  the  method,  the  wind 
effects  on  a  planar  representation  of  a  spray  curtain  are 
exami ned . 


Spray  curtains  are  currently  used  for  the  dispersion 
of  heavy,  sometimes  toxic  or  flammable,  gases.  As  a  step  toward 
evaluating  the  effectiveness  of  the  spray  curtain  for  gas  dis¬ 
persal,  the  flow  field  produced  by  the  interaction  of  the  wind 
with  the  spray  is  predicted. 

4.2.1  Model^l2n£  ^f_a_Plana£  represen;ta;tion 
of_a_S£ray_curta2n 

A  spray  curtain  is  formed  by  distributing  a  series 
of  spray  nozzles  along  a  certain  orientation,  often  forming  a 
straight  line,  and  ideally  situated  normal  to  the  wind  direction 

In  this  orientation,  a  planar  model  of  the  spray 
curtain  is  applied  to  predict  the  global  behaviour  of  the  flow 
field  under  the  influence  of  the  curtain.  Nozzles  of  similar 
type  and  massflow  are  assumed  to  be  equally  distributed  along  a 
straight  line  parallel  to  the  ground.  The  length  of  the  curtain 
should  be  large  compared  with  the  height  of  the  nozzles  so  that 
three  dimensional  effects  are  suppressed. 

In  the  planar  formulation  the  discrete  spacing  of  the 
nozzles  is  represented  as  a  constant  uniform  distribution  of 
injected  liquid  along  the  length  of  a  curtain.  The  distributed 
volume  flow  of  the  curtain  per  unit  length,  q^,  is  calculated  by 
multiplying  the  volume  flow  of  the  individual  nozzles  by  the 
number  of  nozzles  and  dividing  by  the  length  of  the  curtain, 
expressed  by  the  following 
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(4.2.1) 


Here  is  the  number  of  nozzles,  is  the  volume  flow  of 

the  individual  nozzles,  and  is  the  length  of  the  curtain. 


The  spray  in  the  plane  of  symmetry  is  represented  by 
five  to  fifteen  trajectories  of  varying  injection  angle  of 
constant  increment.  The  possibility  of  allowing  the  nozzle  to 
be  tilted  relative  to  the  vertically  downwards  position  by 
degrees  is  accounted  for.  This  is  expressed  below  as 


®i  = 


'^traj'^j 


i  =  l , . .  .  ,N, 


(4,2.2) 


Here  is  the  half  angle  of  the  spray  and  is  the  number 

of  trajectories. 

The  fractional  volume  flow  of  liquid  assigned  to 
each  trajectory  is  a  constant,  being  equally  divided  between 
the  trajectories.  This  may  be  expressed  below 
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(4.2.3) 


Only  constant  droplet  injection  velocities  and  droplet 
sizes  are  considered,  and  may  be  calculated  as  before  using 
equations  (4.1.5)  and  (4.1.2). 

The  particular  curtain  configuration  presented  here 
is  modelled  after  typical  configurations  used  in  the  Buxton 
full-scale  trials  (Ref.  17).  In  these  trials  a  Protectospray 
nozzle  was  used,  but  due  to  the  lack  of  data  on  the  nozzle 
the  following  hypothetical  nozzle  characteristics  are  used. 


Inside  diameter 

Nozzle  design 

Half  angle 

parameter 

of  spray 

^0 

®MAX 

1/3  2/3 

mm 

Kg  /sec 

degrees 

12.7 

.75 

45 

Table  4.2.1  Hypothetical  nozzle  characteristics 

The  volume  flow  through  each  nozzle  is  assumed  to 
obey  the  following  relationship  to  the  delivery  pressure. 

Qw  =  *1334  VaP^  (4.2.4) 

where  aPj^  is  given  in  kP^  and  in  liters  per  second.  Typical 
overpressures  of  127  and  650  kPj  yield  volume  flows  at  1.5  to 
3.4  liters  per  second  through  each  nozzle. 

In  the  Buxton  trials,  curtain  heights  ranged  from  two 
to  three  meters  with  a  nozzle  spacing  of  0.5  to  3  nozzles  per 
meter  covering  distances  at  16  to  32  meters.  These  typical 
values  serve  as  inputs  for  the  calculations  made  here. 

The  computations  are  made  in  a  domain  similar  to 
that  shown  in  figure  4.2.1.  The  nozzle  is  located  approximately 
in  the  middle  of  the  domain. 


The  applied  boundary  conditions  are  also  shown  in 
figure  4.2.1.  At  the  specified  inflow  boundary  a  constant 
horizontal  velocity  and  zero  vertical  velocity  are  applied. 

Optimization  of  the  relaxation  coefficient  for  these 
calculations  yields  a  value  of  1.7. 

The  effective  viscosity  is  taken  as  the  molecular 
viscosity  of  air. 
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4.2.2  Re^u 

The  method  is  applied  to  examine  the  effect  of  wind 
speed,  nozzle  spacing,  and  tilting  the  nozzle  into  or  away  from 
the  wind  on  the  performance  of  a  spray  curtain.  The  effect  of 
curtain  height  and  spray  angle  on  the  curtain  performance  as 
well  as  the  sensitivity  of  the  calculation  to  drop  size  and 
injection  velocity  are  also  examined. 

4. 2. 2.1  Effectof  wind_sgeed 

Examination  of  the  effect  of  wind  speed  on  spray  cur¬ 
tain  performance  may  be  made  by  performing  calculations  on  a 
particular  curtain  geometry  for  increasing  wind  speed. 

The  curtain  examined  here  is  of  2.5  meters  height 
with  the  nozzles  hypothetically  spaced  one  per  meter  and  opera¬ 
ting  at  individual  volume  flow  rates  of  three  liters  per  second. 

Figure  4. 2. 2-4. 2. 5  shows  the  resulting  flow  field 
for  increasing  wind  speeds  from  2  to  4.5  meters  per  second. 

To  be  noted  is  the  decrease  in  the  size  of  the  recirculation 
region  in  front  of  the  spray  with  increasing  wind  speed,  and 
its  complete  disappearance  at  the  higher  wind  speed.  These 
results  demonstrate  well  the  degradation  of  curtain  performance 
with  increasing  wind  speed. 

4. 2. 2. 2  Effect_of  nozzle_spaci ng 

Some  of  the  effects  of  nozzle  spacing  may  be  examined, 
though  three  dimensional  effects,  such  as  interaction  between 
adjacent  sprays  and  flow  between  adjacent  sprays,  are  not 
considered . 

Figures  4.2.6  and  4.2.7  show  the  curtain  with  the 
same  height  and  nozzle  operating  conditions  but  with  a  spacing 
increased  by  a  factor  of  two. 
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Again  the  decrease  in  curtain  performance  is  observed 
with  increasing  wind  speed.  Also  seen  is  the  decrease  in  per¬ 
formance  with  increasing  spacing  when  compared  with  the  previous 
set  of  results.  It  must  be  stated,  however,  that  this  is  true 
as  long  as  the  droplets  from  adjacent  sprays  do  not  interact. 

Figure  4.2.8  and  4.2.9  show  the  curtain  with  the  same 
height  of  2.5  meters  but  with  the  nozzles  spaced  at  one  per 
meter  and  operating  at  a  volume  flow  rate  of  1.5  liters  per 
second.  Comparison  between  these  two  figures  and  the  two  pre¬ 
vious  ones  yields  the  conclusion  that  better  performances  may 
be  obtained  for  the  same  liquid  volume  flow  per  unit  length  by 
spacing  the  nozzles  further  apart  and  operating  each  nozzle  at 
a  higher  massflow.  This  conclusion  is  valid  to  a  limit  where 
the  spacing  between  nozzles  becomes  large  enough  that  the  gas 
slips  through  the  gaps  in  the  curtain.  This  limiting  behaviour 
is  not  predictable  with  this  two-dimensional  model. 

4. 2, 2. 3  |^fect_of  tilting  nozzle 

During  the  Buxton  trials,  it  was  observed  that  tilting 
the  spray  into  the  wind  yielded  better  dispersion  of  a  heavy  gas 
compared  with  the  dispersion  with  the  nozzle  positioned  verti¬ 
cally  downwards  or  tilted  away  from  the  wind. 

Figures  4.2.10-4.2.13  show  the  changes  in  the  flow 
field  as  the  nozzle  is  tilted  increasingly  more  into  the  wind. 

In  the  vertically  downwards  position  there  is  little  modifica¬ 
tion  of  the  wind  field  due  to  the  spray.  As  the  tilt  angle  is 
increased  Into  the  wind  a  stagnation  region  behind  the  spray 
begins  to  appear.  The  flow  upstream  must  then  rise  over  the 
curtain  or  pass  beneath  the  stagnation  region  where  it  seems 
to  be  mixed  somewhat  with  the  air  entrained  by  the  nozzle.  The 
degree  to  which  mixing  between  a  heavy  gas  layer  and  surrounding 
air  might  occur  is  not  immediately  assessable  without  a  disper¬ 
sion  model,  but  it  is  certainly  clear  that  comparatively  stronger 
modification  to  the  flow  field  is  obtainable  by  tilting  the 


nozzle  into  the  wind.  Efforts  toward  coupling  a  dispersion 
model  with  a  flow  field  calculation  are  in  progress  (Ref.  18). 

Figures  4.2.14  and  4.2.15  show  the  nozzle  tilted 
away  from  the  wind.  Here  it  seems  the  flow  is  simply  accelera¬ 
ted  downstream  with  very  little  mixing  occurring.  This  is  cin- 
sistent  with  the  results  of  the  Buxton  trials  where  it  was 
found  that  tilting  the  nozzle  away  from  the  wind  led  to  de¬ 
creased  dispersion  compared  with  the  vertically  downwards  posi¬ 
tion  or  tilted  into  the  wind.  Though  convergence  difficulties 
prevented  obtaining  results  for  higher  volume  flow  rates  per 
unit  length,  indications  of  a  strong  recirculation  downstream 
of  the  nozzles  exist. 

4. 2. 2. 4  Effect_of  decreasi ng_cone  angle 

Figure  4.2.16  shows  the  effect  decreasing  the  cone 
angle  of  the  spray  when  compared  with  figure  4.2.4.  It  appears 
locally  that  the  effect  of  the  spray  is  more  concentrated  for 
the  smaller  cone  angle,  though  globally,  there  is  not  much 
difference  between  the  two  cases. 

4. 2. 2. 5  Effect_of  heightof  spray 

Figure  4.2.17  and  4.2.18  when  compared  with  figure 
4.2.4  allow  the  effect  of  the  height  of  the  spray  to  be  examined 
It  can  be  seen  that  raising  the  height  of  the  spray  extends  its 
region  of  influence,  though  the  increased  spray  height  also 
causes  the  higher  velocities  near  the  nozzle,  which  might  encou¬ 
rage  more  rapid  mixing,  to  be  displaced  further  from  the  ground 
level  (and  possible  further  from  a  heavy  gas  layer). 

The  optimal  curtain  height  it  seems  should  be  large 
enough  so  that  a  heavy  gas  cloud  would  be  contained  by  the  cur¬ 
tain,  but  low  enough  so  that  more  rapid  mixing  between  the  heavy 
gas  and  the  surrounding  air  occurs. 
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4. 2. 2. 6  Sensitivity  of  thecal cul ationto  drop 
sizeandinjectionvelocity 

Figure  4.2.19  and  figure  4.2.20  show  the  flow  fields 
when  droplet  size  is  decreased  by  20%  and  40%  respectively 
compared  to  the  results  shown  in  figure  4.2.4.  Small  variations 
in  the  flow  field  are  noticed.  The  perturbations  introduced  in 
the  main  flow  field  by  the  spray  increase  slightly  with  decreas¬ 
ing  droplet  size. 

Figure  4.2.21  and  4.2.22  depict  the  effect  of  increas¬ 
ing  the  droplet  injection  velocity  by  20%  and  40%  over  the  cor¬ 
responding  injection  velocity  used  to  obtain  the  results  in 
figure  4.2.4.  Seen  here  is  the  increase  in  the  perturbation  of 
the  flow  field  with  increasing  injection  velocity,  with  the 
reappearance  of  the  recirculation  region  occurring  for  the 
higher  injection  velocity. 

4.2.3  Conclusions 


The  model  is  shown  to  give  qualitatively  good  results. 

The  predictions  demonstrate  the  degradation  of  curtain 
performance  with  increasing  wind  speed  and  the  improvements  in 
performance  with  increased  liquid  massflow  through  the  nozzles. 

For  the  same  liquid  volume  flow  per  unit  length,  the 
results  predict  better  performance  by  spacing  the  nozzles  further 
apart  and  operating  each  nozzle  individually  at  a  higher  volume 
flow  rate.  This  result  is  true  to  a  limit  where  the  gas 
begins  to  slip  through  adjacent  sprays,  which  is  not  predictable 
with  this  model . 

Quite  significantly,  the  model  indicates  improved 
curtain  performance  can  be  obtained  by  tilting  the  nozzle  into 
the  wind  direction,  while  tilting  of  the  nozzle  away  from  the 
wind  is  not  advised. 


Results  also  show  that  the  height  of  the  curtain  can 
be  optimized  so  that  a  heavy  gas  cloud  does  not  roll  over  the 
curtain,  while  ensuring  that  entrained  air  has  sufficient  energy 
to  mix  with  the  heavy  gas  cloud. 

Extension  of  the  method  to  include  a  dispersion  model 
is  recommended  and  is  currently  proceeding. 

4 . 3  Calculation  of  the  flow  inside  a  cooling  tower 

The  method  is  applied  to  calculate  the  velocity  and 
pressure  fields  inside  and  near  the  inlet  of  a  cooling  tower. 

Of  particular  interest  is  the  effect  of  the  droplets  on  the 
velocity  and  pressure  fields. 

4.3.1  Mo d e  1_1  j[n£  £ f_t he_coo  l_i  ng_t owe r 

The  domain  of  interest  is  primarily  the  lower  portion 
of  the  cooling  tower  below  the  packing  which  contains  the  water 
droplets.  A  computational  domain  (shown  in  Fig.  4.3.1)  is  con¬ 
structed  to  allow  the  the  investigation  to  be  focussed  in  this 
region.  The  computational  domain  is  equally  divided  between  the 
inside  and  outside  of  the  tower  so  that  the  flow  entering  the 
tower  would  not  be  greatly  affected  by  the  flow  at  the  outer 
boundaries . 


The  shell  of  the  tower  is  represented  by  a  hollow 
cylinder  of  zero  thickness  and  of  radius  located  at  a 

distance  above  the  ground. 


The  packing,  of  thickness  Lp,  is  located  a  distance 
of  above  the  ground  and  extended  across  the  entire  radius 

fhA  Tho  narlfinn  i  c  rnnHelloH  iicinn  an  annnnarh  h  v/ 
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of  the  tower.  The  packing  is  modelled  using  an  approach  by 
Buffat  et  al .  (Ref.  8)  which  treats  the  packing  as  a  porous 
media  that  causes  a  loss  of  momentum  to  the  gas  passing  through 
it.  The  momentum  loss  per  unit  volume  may  be  expressed  in  the 
following  manner  : 


40 


K 

~  (4.3.1) 

K2 

F^  - - lV^(r,z)  lV^(r,z)  (4.3.2) 

These  expressions  when  adapted  to  the  SMAC  cell  structure  may 
be  written  as  follows  ; 

Ki 

F.  =  —  I  V  (4.3.3) 

i  +  l/2,j  Lp  ^i  +  l/2,j  '^i+l/2,j 

K 

=  -  —  |V  I  V  (4.3.4) 

2iJ  +  l/2  L  ^i.j  +  1/2  ^i,j  +  l/2 


The  velocity  V  ,  not  explicitly  a va i 1 abl e  , '  i  s 

obtained  by  simple  averaging  of  the  four  surrounding  axial 
velocity  components. 

The  values  of  the  coefficients  Ki,K2  are  both  taken 
as  equal  to  6,  as  suggested  in  reference  8. 

The  boundary  conditions  applied  at  the  edges  of  compu¬ 
tational  domain  and  along  both  sides  of  the  tower  shell  are  dis¬ 
played  in  figure  4.3.1.  Along  the  specified  inflow  boundary  a 
constant  normal  velocity,  proportional  to  the  volume  flow  of  air 
through  the  tower,  is  applied.  The  tangential  velocity  along 
the  boundary  is  set  equal  to  that  just  upstream.  The  details  of 
the  boundary  conditions  at  the  other  boundaries  have  been  dealt 
with  in  previous  sections  and  are  not  discussed  here. 

Computations  are  made  on  a  20  x  20  mesh  with  the  drop¬ 
let  phase  represented  by  ten  trajectories  spaced  at  equal  inter¬ 
vals  across  the  packing.  The  volume  flow  of  liquid  assigned  to 
each  trajectory  is  such  that  the  total  volume  flow  per  unit  area 
leaving  the  packing  is  a  constant  (in  a  discrete  sense). 
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This  may  be  expressed  as  follows 


H. 
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nl et 


tower 


N 


tra  j 


(i-1) 


(4.3.5) 
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tra  j 


where 


and  Rj  are  the  vertical  and  horizontal  locations  of 


droplet  injection,  f»  is  the  fractional  volume  flow  of  liquid 


assigned  to  the  i^*^  trajectory  of  the  total 


Qw-j 


N.  .  trajectories. 

t  rd  J 


The  droplet,  of  sizes  varying  from  3  to  6  millimeters, 
are  assumed  to  drip  from  the  lower  edge  of  the  packing  with  an 
initial  vertical  trajectory  and  a  prescribed  initial  vertical 
velocity. 


Uq.  =  0.2  m/s  (4.3.6) 

e.  =  0 


The  effective  viscosity  of  the  fluid  is  taken  as 
1  m  /s  to  account  for  the  turbulence  induced  by  the  droplets 
(except  where  noted). 

The  relaxation  coefficient  for  the  pressure-velocity 
modification  is  chosen  as  1.4,  as  used  for  computation  of  axisym- 
metric  spray,  though  no  optimization  of  this  parameter  is  made 
for  the  case  of  a  cooling  tower. 
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4.3.2 

The  method  is  applied  to  a  cooling  tower  with  dimen¬ 
sions  shown  in  Table  4.3.1.  Calculation  of  the  velocity  and 
pressure  fields  is  made  for  the  cases  of  no  droplets  and  drop¬ 
lets  of  size  3,  4  and  6  millimeter  diameter. 


Tower  radius 

p 

tower 

60 

m 

Inlet  height 

^  i  n  1  e  t 

10 

m 

Packing  thickness 

s 

2 

m 

Volume  flow  of  air 
through  tower 

Qa 

22619 

m^/s 

Volume  flow  of  water 

31.4 

m  Vs 

Table  4.3.1  Tower  dimensions  and  test  conditions 

A  volume  flow  of  air  inside  the  tower  equal  to 
22619  m^/s  is  used  yielding  a  mean  velocity  through  the  packing 
of  2  m/s.  For  the  cases  considered  with  droplets,  a  volume  flow 
of  water  equal  to  31.4  m^/s  gives  a  constant  volume  flow  per 
unit  are  of  10  m^/hr-m^. 

4.3.2. 1  Veloci ty_f i eld 

The  computed  velocity  fields  for  the  above  conditions 
with  no  droplets,  and  droplets  of  3,  4  and  6  millemeters  are 
shown  in  figures  4.3.2,  4.3.4,  4.3.6  and  4.3.8,  respectively. 

The  effect  of  the  packing  on  the  flow  field  is  clear. 
The  radial  component  of  velocity  is  suppressed  significarttly  as 
gas  passes  through  the  packing,  while  the  axial  component  is 
caused  to  be  fairly  uniform  (±  15%  variation)  across  the  radius 
of  packing  except  very  near  the  tower  shell  where  the  velocity 
goes  to  zero. 


The  effect  of  the  droplets  on  the  velocity  field  is 
to  cause  a  retardation  of  penetration  of  the  gas  towards  the 
axis  of  symmetry  though  this  effect  is  somewhat  dominated  by 
the  opposing  effect  of  the  packing  which  tends  to  cause  the 
axial  velocity  to  be  radially  uniform  across  the  packing.  The 
effect  is  somewhat  more  clearly  shown  in  figures  4.3.10  and 
4.3.11. 

Additionally  the  effect  of  the  droplets  is  to  cause 
the  profiles  of  the  radial  component  of  velocity  below  the  pack¬ 
ing  to  become  slightly  more  uniform,  thus  decreasing  the  size 
of  the  stagnation  region  in  the  lower  corner  of  the  domain  near 
the  axis  of  symmetry. 

4. 3. 2. 2  Pressure_fie1d 

The  corresponding  pressure  fields  for  the  above  compu¬ 
tations  are  shown  in  figures  4 . 3 i 3 , '4 .3 . 5 ,  4.3.7  and  4.3.9. 

Here  non  dimensional  static  pressure  isobars  are  displayed  in 
increments  of  1.2.  Pressure  has  been  nondimensional  ized  by  the 
gas  density  times  the  square  of  the  mean  axial  velocity  through 
the  packing  (2  m/s) . 

In  these  figures  the  pressure  losses  caused  by  the 
packing  and  the  inlet  are  clearly  shown  by  the  crowding  of  iso¬ 
bars  in  these  regions. 

The  pressure  loss  due  to  the  droplets  is  shown  to  be 
quite  significant  and  increases  with  decreasing  droplet  diameter 
The  physical  validity  of  this  trend  can  be  demonstrated  by  the 
following  argument.  The  pressure  drop  due  to  the  droplets 
is  proportional  to  the  product  of  the  number  of  droplets  and 
the  average  droplet  drag.  The  number  of  droplets  inside  the 
tower  can  be  expressed  in  terms  of  the  volume  flow  of  water,  the 
time  for  the  droplet  to  fall  from  the  packing  to  the  ground 
and  the  droplet  diameter,  and  the  drag  in  terms  of  the  drag 
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coefficient,  droplet  velocity  relative  to  the  air,  and  droplet 
diameter.  Thus  the  pressure  drop  due  to  droplets  may  be  ex¬ 
pressed 


IT 
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Noting  that  the  forces  acting  on  the  droplet  are  dominated  by 
the  gravitational  force,  the  time  for  the  droplet  to  reach  the 
ground  and  the  relative  velocity  of  the  droplet  may  be  assumed 
to  be  mild  functions  of  droplet  diameter  (decreasing  and  in- 
creasily  respectively). 

Further  neglecting  the  small  variations  in  the  drag 
coefficient  with  droplet  diameter,  the  pressure  loss  is  seen  to 
vary  with  the  inverse  of  the  droplet  diameter  for  constant 
volume  flow  of  water. 

Calculation  of  the  additional  total  pressure  loss 
caused  by  the  droplets,  shown  in  table  4.3.2,  is  conventionally 
obtained  by  subtracting  the  total  pressure  loss  with  droplets 
from  the  total  pressure  loss  in  the  absence  of  droplets.  The 
total  pressure  loss  in  each  case  is  obtained  by  subtracting 
the  average  total  pressure  above  the  packing,  from  the 

average  total  pressure  outside  the  tower  far  from  the  inlet 
(see  figure  4.3.1). 

Further  assuming  the  pressure  loss  due  to  the  packing 
to  be  equal  to  the  coefficient  K2  times  the  average  axial  velo¬ 
city  through  the  packing  squared,  allows  the  approximate  pressure 
loss  due  to  the  inlet  to  be  extracted.  These  are  also  displayed 
in  table  4.3.2. 

The  calculated  pressure  loss  due  to  the  droplets  is 
reasonable  though  experimental  data  seem  to  indicate  that  these 
values  slightly  underestimate  the  actual  pressure  loss.  Further 
experimental  verification  is  warranted. 


1 


Overall  total  pressure 

losses 

(Pt 

■2-Pt 

1) 

No  droplets 

7.6 

mm 

H2O 

D  =  3  mm 

10.7 

mm 

H2O 

I 

D  =  4  mm 

10.1 

mm 

H2O 

D  =  6  mm 

9.4 

mm 

H  0 

1 

Individual  total 

pressure  losses 

1 

1 

Due  to  inlet 

4.6 

mm 

H20 ! 

Due  to  packing 

-^3.0 

mm 

H20 

Due  to  droplets 

D  = 

3  mm 

3.1 

mm 

H20 

D  = 

4  mm 

2.5 

mm 

H20 

D  = 

6  mm 

1.8 

mm 

H20 1 

Table  4.3.2  Total  pressure  losses  in  domain 
of  interest 


Overall  total  pressure  losses  (Pt2"*’ti) 

No  droplets  5.8  mm  H2O 

D  =  4  mm  8.4  mm  H2O 

Individual  pressure  losses 
Due  to  inlet  2.8mmH20 

Due  to  packing  '>-3.0  mm  H2O 

Due  to  droplets  D  =  4  mm  2.65  mm  H2O 


Table  4.3.3  Total  pressure  losses  in  domain 
of  interest  =  0.1  m^/s 
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The  total  pressure  loss  due  to  the  inlet,  however,  is 
grossly  overestimated.  This  is  most  probably  due  to  the  high 
effective  viscosity  chosen  to  represent  the  turbulence  induced 
by  the  droplets,  which  may  not  model  the  effective  viscosity 
outside  the  inlet  where  there  are  no  droplets. 

This  is  substantiated  by  the  calculations  shown  in 
figures  4.3.12,  4.3.13,  4.3.14  and  4.3.15  in  which  an  effective 
viscosity  of  0.1  m^/s  was  used.  The  overall  and  individual 
pressure  losses  shown  in  table  4.3.3  indicate  a  large  reduction 
in  the  pressure  drop  due  to  the  inlet,  which  demonstrates  the 
sensitivity  of  this  component  of  the  pressure  loss  to  the  effec¬ 
tive  viscosity.  This  indicates  that  a  variable  effective  visco¬ 
sity  may  be  required  to  properly  simulate  the  different  nature 
of  the  flow  inside  and  outside  the  cooling  tower. 

Figure  4.3.16  41so  indicates  differences  in  the  velo¬ 
city  field.  Compared  here  are  the  velocity  profiles  at  the  in¬ 
let  and  after  the  packing  for  the  case  of  4  mm  droplets  and  ef¬ 
fective  viscosities  of  1.0  and  0.1  m^/s.  The  velocity  profiles 
become  flatter  as  the  effective  viscosity  is  decreased. 

Because  the  velocity  and  pressure  fields  are  sensitive 
to  the  effective  viscosity,  further  work  with  regard  to  a  tur¬ 
bulence  modelling  is  recommended. 

4.3.3  £o£cJ[u£ion^ 

The  results  demonstrate  the  applicability  of  the  method 
for  the  computation  of  flow  through  a  cooling  tower.  Results  are 
reasonable  with  the  exception  of  the  pressure  drop  due  to  the  in¬ 
let.  Improvement  is  expected  with  further  work  in  the  modelling 
of  the  turbulence. 

Extension  of  the  method  to  account  for  the  true  tower 
geometry  is  possible  as  the  SMAC  method  is  applicable  to  problems 
with  curved  boundaries  (Ref.  5). 
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5.  CONCLUSIONS 

The  results  demonstrate  the  successful  implementation 
of  the  model  with  good  qualitative  agreement  being  seen  through¬ 
out.  Comparison  with  experimental  results  indicated  a  slight 
underprediction  of  the  effect  of  the  droplets  though  results 
are  quite  reasonable  and  demonstrate  the  usefulness  of  the 
model . 

Application  of  the  model  to  axisymmetric  sprays 
demonstrates  the  limited  effect  of  the  ceiling  on  spray  perfor¬ 
mance,  and  the  lack  of  influence  of  the  pressure  term  on  the 
flow  field  prediction  validating  some  assumptions  made  in  other 
investigations. 

The  calculations  for  the  planar  representation  of  a 
spray  curtain  demonstrate  the  degrading  effect  of  wind  on  its 
performance,  but  indicate  that  this  may  be  overcome  by  tilting 
the  nozzle  into  the  wind. 

The  modelling  of  the  flow  in  a  cooling  tower  demon¬ 
strates  the  reasonable  prediction  of  the  effects  of  the  droplets 
on  both  the  velocity  and  pressure  fields. 

Further  development  of  the  model  is  recommended,  par¬ 
ticularly  in  the  areas  of  turbulence,  dispersion  and  heat 
transfer  model  ling. 
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DOMAIN  BOUNDARY 


FIG.3.1-SMAC  CELL  STRUCTURE  FOR  THE  GAS  PHASE  COMPUTATION. 


- ►  CORRESPONDS  TO  A  GAS  VELOCITY  =  2-20  M/S 

SZl  NOZZLE  30.  DEGREE  HALF  ANGLE 

VOLUME  FLOW  OF  WATER  =  6-6  LITERS  PER  MINUTE 

DROPLET  injection  VELOCITY  =  7.1  M/S 

MEAN  DROP  SIZE  =  -83  MM  HEIGHT  OF  NOZZLE  =  2-0  M 


FIGURE  4-1-2  -PLOW  FIELD  DUE  TO  A  SPRAY 


- ►  CORRESPONDS  TO  A  GAS  VELOCITY  =  5-00  M/S 

SZl  NOZZLE  30-  DEGREE  HALF  ANGLE 

VOLUME  FLOW  OF  WATER  =  13-3  LITERS  PEK  MINUTE 

DROPLET  INJECTION  VELOCITY  =  14-6  M/S 

MEAN  DROP  SIZE  =  • 48  MM  HEIGHT  OF  NOZZLE  =  2-0  M 

FIGURE  4-1-3  -FLOW  FIELD  DUE  TO  A  SPRAY 


- -  CORRESPONDS  TO  A  GAS  VELOCITY  =  8-00  .M/S 

SZl  NOZZLE  30.  DEGREE  HALF  ANGLE 

VOLUME  FLOW  OF  WATER  =21-1  LITERS  PER  MINUTE 

DRO.PLET  INJECnON  VELOCITY  =  23-2  M/S 

MEAN  DROP  SIZE  -  .32  MM  HEIGHT  OF  NOZZLE  =  2-0  M 

FIGURE  4-1.4  -FLOW  FIELD  DUE  TO  A  SPRAY 


- -  CORRESPONDS  TO  A  GAS  VELOCITY  =  4.10  M/S 

SZl  NOZZLE  30.  DEGREE  HALF  ANGLE 

VOLUME  FLOW  OF  WATER  =  6.6  LITERS  PER  MINUTE 

DROPLET  INJECTION  VELOCITY  =  7-2  M/S 

MEAN  DROP  SIZE  =  .83  MM  HEIGHT  OF  NOZZLE  =  0-4  M 


FIGURE  4-1-5  -FLOW  FIELD  DUE  TO  A  SPRAY 


- *  CORRESPONDS  TO  A  GAS  VELOCITY  =  9.50  h/S 

SZl  NOZZLE  30.  DEGREE  HALF  ANGLE 

VOLUME  flow  of  WATER  =  13.3  LITERS  PER  MINUTE 

DROPLET  INJECTION  VELOCITY  =  M.6  M/S 

MEAN  DROP  SIZE  =  •  ■48  MM  HEIGHT  OF  NOZZLE  =  0-4  M 

FIGURE  4-1-6  -FLOW  FIELD  DUE  TO  A  SPRAY 


- ►  CORRESPONDS  TO  A  GAS  VELOCITV  -  1 7 ■ 00  f 

SZ1  NOZZLE  30.  DEGREE  HALF  ANGLE 

VOLUME  FLOW  OF  WATER  -  21-1  LITERS  PER  MINUTE 

DROPLET  INJECTION  VELOCITY  =■•  23-2  M/S 

MEAN  DROP  SIZE  -  -32  MM  HEIGHT  OF  NOZZLE  =  2-0  M 


FIGURE  4-1-7  -FLOW  FIELD  DUE  TO  A  SPRAY 


MEASUREMENT  SPRAYS 


-PERFORMANCES  DATA  FOR  SPRAYS 


S22  N022LE  30.  DEGREE  HALF  ANGLE 


ERSUS  INVERSE  SPRAY  DIAMETER  SQUARED 


RADIAL  COORDINATE  (CM) 

FIGURE  4-1-12  -PREDI  CTED  AND  MEASURED  VELDCITIES 
INSIDE  THE  SPRAY 


SZ2  NOZZLE 
C  =  44-5  LT/MIN 
Z=  40.  CM 


F1 CURE ;  4. 1  .16  -CQMPARl SON  OF  RESULTS  WITH  AND  WITHOUT  PRESSURE  TERMS 


UNSPECIFIED 


N 


AdVONnOa  ANDIJNI  a3UID3dS 


FIG.  A.  2.1 -COMPUTATIONAL  DOMAIN  FOR  PLANAR  SPRAY  CURTAIN 
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